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Abstract 

Recently, Shi and Sun proposed Point Integral method (PIM) to discretize Laplace-Beltrami operator 
on point cloud [161 il9j . In PIM, Neumann boundary is nature, but Dirichlet boundary needs some special 
treatment. In our previous work, we use Robin boundary to approximate Dirichlet boundary. In this 
paper, we introduce another approach to deal with the Dirichlet boundary condition in point integral 
method using the volume constraint proposed by Du et.al. [?]■ 


1 Introduction 


Partial differential equations on manifold appear in a wide range of applications such as material science 
mm, fluid flow mm, biology and biophysics Emnnnna and machine learning and data analysis 01 El- 
Due to the complicate geometrical structure of the manifold, it is very chanlleging to solve PDEs on manifold. 
In recent years, it attracts more and more attentions to develop efficient numerical method to solve PDEs 
on manifold. In case of that the manifold is a 2D surface embedding in R 3 , many methods were proposed 
include level set methods mm, surface finite elements [8], finite volume methods m, diffuse interface 
methods m and local mesh methods [T3] , 

In this paper, we focus on following Poisson equation with Dirichlet boundary condition 

f -A m u(x) = /(x), xeM , , 

\ u(x) = 0, X € dM [ ' 1 

where M is a smooth manifold isometrically embedded in R d with the standard Euclidean metric and dM 
is the boundary. A m is the Laplace-Beltrami operator on manifold M. Let g be the Riemannian metric 
tensor of M. Given a local coordinate system (x 1 ,# 2 , ■ ■ ■ ,x k ), the metric tensor g can be represented by a 
matrix [gij]kxk, 

d d . . 1 

9ij < dx i 'dxi > ' M 


Let [g lJ ]kxk is the inverse matrix of [gij\kxk, then it is well known that the Laplace-Beltrami operator is 


A m = 


d 


d 




In this paper, the metric tensor g is assumed to be inherited from the ambient space R d , that is, M 
isometrically embedded in R d with the standard Euclidean metric. If M is an open set in R d , then Am 
becomes standard Laplace operator, i.e., A m — 5b=i jp?- 
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In our previous papers, mm, Point Integral method was developed to solve Poisson equation in point 
cloud. The main observation of the Point Integral method is that the solution of the Poisson equation can 
be approximated by an integral equation, 

7 / Rt(x,y)(u(x)-u(y))dy-2 [ #t(x,y)|^(y)d/r y = f R t (x, y)/(y)dy (1.2) 

£ J M JdM aYl J M 

where n is the out normal of M. at dM. The kernel functions 


#t(x, y) 


i pf ll x -y|H 

(47rt) fe / 2 \ 4 1 ) ’ 


#t(x, y) 


i ^ iN-yin 

(47rt) fc / 2 \ 4 1 ) 


(1.3) 


and R(r) = f+°° R(s)ds. t is a parameter, which is determined by the desensity of the point cloud in the 
real computations. 

The kernel function R(r) : R + —»■ R+ is assumed to be C 2 smooth and satisfies some mild conditions (see 
Section mi) . 

The integral approximation (11.21) is natural to solve the Poisson equation with Neumann boundary 
condition. To enforce the Dirichlet boundary condition, in our previous work os na, we used Robin 
boundary condition to approximate the Dirichlet boundary condition. More specifically, we solve following 
problem instead of EH with 0 < /J < 1, 


-A M u(x) = /(x), x€M, 

u(x) + j3§£ = 0, xedM. 

Using El, we have an integral equation to approximate the above Robin problem, 


/ #t(x,y)(u(x)-u(y))dy + - / R t (x,y)u(y)dfiy = R t (x, y)/(y)dy. 
>M P JdM JM 


(1.4) 


(1.5) 


We can prove that this approach converge to the original Dirichlet problem |T9]. In the real computations, 
small fj may give some trouble. The overcome this problem, we also introduced an itegrative method to 
enforce the Dirichlet boundary condition based on the Augmented Lagrangian Multiplier (ALM) method. 
However, we can not prove the convergence of this iterative method, although it always converges in the 
numerical tests. 

Recently, Du et.al. 17] proposed volume constraint to deal with the boundary condition in the nonlocal 
diffusion problem. They found that in the nonlocal diffusion problem, since the operator is nonlocal, only 
enforce the boundary condition on the boundary is not enough, we have to extend the boundary condition 
to a small region close to the boundary. Borrowing this idea, in nonlocal diffusion problem to handle the 
Dirichlet boundary. This idea gives us following integral equation with volume constraint: 


Rt(x,y){u(x) — u(y))dy = / R t (x, y)/(y)dy, x £ M' t 




u(x) = 0, 


JM 


( 1 . 6 ) 


x € V 4 


Here, M! t and V t are subsets of M. which are defined as 

M' t = jx e A4 : B ^x, 2\ft^ D dM. = 0 j , V t = M\M' t . (1.7) 

The thickness of V* is 2 \Jt which implies that |V*| = 0{\fi). The relation of A4, dM, M' t and Vt are sketched 
in Fig. Q] 

The main advantage of the integral equation (11.61) is that there is not any differential operator in the 
integral equation. Then it is easy to discretized on point cloud. Assume we are given a set of sample points 
P = {p, : p; £ M, i = 1, ■ ■ • , n} sampling the submanifold M and one vector V = (Vi, ■ ■ • , 14) 4 where V) 
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M = M' t U Vt 



is the volume weight of p; in Ad. In addition, we assume that the point set P is a good sample of manifold 
Ad in the sense that the integral on A4 can be well approximated by the summation over P, see Section [lj] 
Then, (11.611 can be easily discretized to get following linear system 

R t(Pi>Pj)( u i- u j)V j = p t eA4' t , 

1 p Pj€ M (1-8) 

Ui = 0, Pi G Vi. 

This is the discretization of the Poisson equation m given by Point Integral Method with volume constraint 
on point cloud. 

Similarly, the eigenvalue problem 


—A m u ( x ) = Aw, x G Ad 
u(x) =0, x e dAi 


can be approximated by an integral eigenvalue problem 

1 


t / #t(x,y)(u(x)-u(y))dy = A / P t (x, y)u(y)dy, x G M' t 

1 J M JM 


u(x) = 0, 

And corresponding discretization is given as following 


x € Vi 


Rt(Pi,Pj){ui - Uj)Vj = A ^2 Pt(Pi, Pj)wjV j, P i&M't, 


P j&M 


PjGM 


Ui = 0 , 


p* e v t . 


(1.9) 


( 1 . 10 ) 


( 1 . 11 ) 


1.1 Assumptions and main results 

One of the main contribution of this paper is that, under some assumptions, we prove that the solution of 
the discrete system (11.81) converges to the solution of the Poisson equation (11.11) and the spectra of the eigen 
problem (11.111) converge to the spectra of the Laplace-Beltrami operator with Dirichlct boundary (11.91) . 

The assumptions we used are listed as following. 

Assumption 1.1. • Assumptions on the manifold: A4,dA4 are both compact and C°° smooth. 

• Assumptions on the sample points (P, V): (P, V) is h-integrable approximation of Ad, i.e. 
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For any function f £ C 1 (M), there is a constant C independent of h and f so that 


IM 


E 

PiGM 


/(y)dy- V f(Pi)Vi 


< Ch\supp(f)\\\f\\ c i {M) . 


• Assumptions on the kernel function R(r): 

(a) R £ C' 2 (M+); 

(b) R(r) > 0 and R(r) = 0 for Vr > 1; 

(c) 3i5o > 0 so that R(r) > So for 0 < r < 

These assumptions are default in this paper and they are omitted in the statement of the theoretical 
results. And in the analysis, we always assume that t and h/y/t are small enough. Here, ’’small enough” 
means that they are less than a generic constant which only depends on A4. 

Under above assumptions, we have two theorems regarding the convergence of the Poisson equation and 
corresponding eigenvaule problem. 

Theorem 1.1. Let u(x) be solution of (ED and u = [mi,-- - ,u n ] 4 be solution of (11.81) and f £ C 1 (A4) in 
both problems. There exists C > 0 only depends on M and dM, such that 

\\ u - < C [t 1/A + II/He 1 (At) 

where 

».»(*) = ! ,>,v, + i Y. 

’■ K ’ I ’ V ' \pjeAt PjGM 

{ 0, x £ V t . 

and w t ,h(x.) = E p jeM R t ( x > Pj) V j- 

Theorem 1.2. Let Xi be the ith largest eigenvalue of eigenvalue problem ED- And let A*’ ,! be the ith largest 
eigenvalue of discrete eigenvalue problem (11. IIP , then there exists a constant C such that 

\Xl’ h - Ai| < CXj [t 1/A + , 

and there exist another constant C such that, for any <f> £ E(Xi,T)X and X = H l (AA' t ), 

||/» - T t , h )<j>\\ HH M' t ) < C {t 1/A + • 

where o\' h = {A*' ,1 £ <r(T tt h) ■ j £ h} and Ii = {j £ N : A j = Ai} ; E(X,T) is the Riesz spectral projection 
associated with X. 


x £ M' t , 


( 1 . 12 ) 


2 Stability analysis 


To prove the convergence, we need some stability results which are listed in this section. The first lemma is 
about the coercivity of the integral operator and the proof can be found in [El- 

Lemma 2.1. For any function u £ L 2 (A4), there exists a constant C > 0 only depends on JXi, such that 




Rt(x, y)(u(x) 


u(y)) 2 dxdy > C I |Vu| 2 dx, 

Jm 
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where 


w ( x ) = —UT / R t{^y)u{y)dy, 
wt(x) J M 

and w t (x) = R t (x, y) dy. 

Jm 

Next corollary directly follows from Lemma [7TT1 

Corollary 2.1. For any function u £ L 2 (M' t ), there exists a constant C > 0 only depneds on Ai, such that 

7 / / ^t(x,y)(u(x) - n(y)) 2 dxdy + j f u 2 (x)([ i? t (x,y)dy^) dx > C [ |Vu| 2 clx, 

1 jM' J M’ 1 JM' \Jvt / jM’ t 


where 


u ( x ) = —/ - R *( x >y)' u (y) d y> 

Wi(x) 7^; 


and wj t (x) = R t (x, y) dy. 

JM 

Proof. Let 


!(x) = 


u(x), xeM' t , 
0, x £ V t . 


Using Lemma T2. 11 


[ |Vn| 2 dx < [ |Vu| 2 dx 

Jm' Jm 


< 


IM 

C 


Rt (x, y)(u(x) - u(y)) 2 dxdy 


IM JM 


7 -f [ -Rt(x,y)(n(x) - n(y)) 2 dxdy + / u 2 (x) f [ i? t (x,y)dy 

t JM' JM' t JM' \Jv t 


dx. 


□ 


Using Lemma EUl we can also get following lemma regarding the stability in L 2 (AA). 

Lemma 2.2. For any function u £ L2(A4) with it(x) = 0 in Vt, there exists a constant C > 0 independent 
on t 

7 f [ #t(x,y)(u(x) — n(y)) 2 dxdy > C\\u\\l 2(M) , 

1 jm Jm 

as long as t small enough. 

Proof. Let 

v(x) = — TA R t(x-,y)u{y)dy. 

wt(x) Jm 

Since u(x) =0, x £ V t , we have 


u(x) = 0, x £ dM. 

By Lemma T2. II and the Poincare inequality, there exists a constant C > 0, such that 

/ |u(x)| 2 dx< f | Vu(x)| 2 dx < y [ f R t (x,y) (it(x) - u(y)) 2 d^ x d^ 
Jm Jm 1 Jm Jm 
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Let S = 75 —-. If u is smooth and close to its smoothed version v, in particular, 

■^L/Irnax 1 W-~:~ 


max 1 co mm 


f w 2 (x)d/r x > 5 2 f u 2 (x)dp 

Jm Jm 


( 2 . 1 ) 


then the proof is completed. 

Now consider the case where EH> does not hold. Note that we now have 

Wu-vWl^m) > ~ IMIl 2 (.M) > 

> = — II»IIm«>- 

0 Wmin 


Then we have 


Ct 

t 


R 


'M JM 


*-y r 


At 


(m(x) - u(y)) 2 dfi x dfi y 


= ^ / «(x) /' R 

IM 


t 


/ M 


l x ~y|- 

At 


(w(x) - w(y))d/i y d/i x 


'M 


'M 


M 2 (x)w t (x)d/i x — / u(x)u(x)u; t (x)d/j, 

J M 

(u(x) - v(x)) 2 wt(x)d[i x + / (u(x) - v(x))v(x)wt (x)d/i 3 

JM 


> 


> 


> 


7/ (m(x) - w(x)) 2 w; t (x)d/i x - 7 f f u 2 (x)«; t (x)d/i x ^ ( f (w(x) - w(x)) 2 wj t (x)d/i : 

£ Jm £ \Jjm / Vx 


1/2 


2 lCn 


[ (w( x ) - w(x)) 2 d^ x - 




2 (x)d/i, 


1/2 


'A4 


(m(x) - v(x.)) 2 dp. 


1/2 


t 


) M 


(m(x) - r(x)) 2 d/r x > -^(1 - <5) 2 / M 2 (x)d/r x . 


/x 


This completes the proof for the theorem. □ 

Corollary 2.2. For any function u G L 2 {M' t ), there exists a constant C > 0 independent on t, such that 


I I #t(x,y)(u(x) - u(y)) 2 dxdy + [ w 2 (x) f j R t (x, y)dy) dx > C\\u\\ 2 L , MI) , 

JM' t JM' t JM' t \Jvt J 


as long as t small enough. 
Proof. Consider 


i(x) = 


u(x), x G M' t , 
0, x G V t . 


and apply Lemma 12.21 

Now, we can prove one important theorem. 

Theorem 2.1. Let u(x) G L 2 (A4) be solution of following integral equation 

1 

t 


□ 


■ / -Rt(x,y)(u(x) - ii(y))dy = r(x) 

JM 

w(x) = 0, 

There exists C > 0 only depends on M. and dM , such that 

INI H l (M' t ) < C\\ r \\ L 2 (M' t ) + Ct\\Vr\\ L 2( M ^ 


xG M' t 
x G V t 


( 2 . 2 ) 
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Proof. First of all, we have 


7 / u(x) / #t(x,y)(u(x)-u(y))dydx 
1 ./ai; J M 

7 / m(x) [ -R t (x,y)(u(x)-w(y))dydx+i / u(x) f i?t(x, y)(u(x) - u(y))dydx 

r -/at; -/a-i; r aa/i; Jv t 

[ [ i?t(x,y)(w(x) — w(y)) 2 dxdy+ j [ w 2 (x) [ i? t (x,y)dydx. 

JA/C J Mi 1 J Mi JVt 


1 

2 t 


/A4J ./At; <- -/A1. 

Now we can get L 2 estimate of u. Using Corollarv l2.21 we have 


Ml2,A4; — 


< 


7 - / [ i?t(x,y)(-u(x)-u(y)) 2 dxdy+y / |w(x )| 2 f f i? t (x,y)dy 

1 JM' t JM' t 1 JM' t \JV t 

u(x) / R t (x,y)(u(x) - u(y))dyd: 

1 J M 


dx 


J M' t J M 

< CHuHa^'lIrlla^ 


This gives that 


IMIl 2 (a/i;) < C|MIl 2 (aij)- 


(2.3) 


Next, we turn to estimate the L 2 norm of Vet in M.' t . Using the integral equation (12.21) . u has following 
expression 


u t [x) = 


1 


wt(x) J M f 

Then ||N7-u,t||| M , can be bounded as following 

1 


i? t (x,y)M t (y)dy + 


t 


w t (x) 


(x), xeVt. 


(2.4) 


IVutU 2 ^, < C 


Wt(x) 


I Mi 


i?t(x,y)ut(y)dy 


+ Cf 2 


2 ,M' t 


r ( x ) 

w t (x) 


2.A4J 


(2.5) 


Corollary 12.11 gives a bound the first term of (12.51) . 


i?t(x,y)u t (y)dy 


2M' t 


< 


7 - [ [ Rt{x,y){u t (x) -w t (y)) 2 dydx+ ^ f \u t (x)\ 2 ( f R t {x, y)dy 

1 J Mi J Mi 1 J Mi \JVt 


( 2 . 6 ) 


dx. 


The second terms of (12.51) can be bounded by direct calculation. 


r(x) 

w t (x) 


C C 


2,M't 


Vr(x) 


w t (x) 


+ C 


2,M't 


r(x)VM f (x) 


(w t (x )) 2 


(2.7) 


2 ,At; 


C 


< CJVr(x)|| 2] _ M , + _ ||r(x )|| 2 _ M , . 

Now we have the bound of || Vitfl^A'i; by combining (12.51) . (12.61) . and (12.71) 

C /' 


Rt(x,y)(u t (x) - u t { y)) dxdy 


’Mi J Mi 


+ J Jm' M x )| 2 M x 1 y)dy^j dx + Ct 2 ||Vr(x )|| 2 _ M , +Ct ||r(> 


( 2 . 8 ) 


‘2,M' t 
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Then the bound of ||Vut|| 2 .xj can be obtained also from (12.81) 

\\Vu t \\l M , 

- T I f - ut(y)) 2dxd y + Ct \\r(x)\\ 2 2M , 

1 JM' t JM' t 

l M ' (X E ^ X,y ^ dy ) dx + Ct2 H Vr ( X )ll2,7KJ 


Then we have 


C 

T 


< 


— [ Wt(x) [ i?t(x,y)(u t (x) - M t (y))dydx 
1 J m ; J m 


+Ct 2 ||Vr(x)||2 j _ M , + Ct \\r{x)\\l M , t 

< \\uth,M't\\ r h,M' t + ct 2 ll Vr (x)ll2,At; + Ct lk( x )ll2,M; 

< C'||r||| 2( ^ ) + Ct 2 ||Vr(x)||^,. 


II VMtlb.THj < C\\r\\ L 2 (M , t) + Ct ||Vr(x )|| 2 _ M , . 
The proof is completed by putting (12.31) and (12.91) together. 


(2.9) 

□ 


3 Convergence analysis 

The main purpose of this section is to prove that the solution of (11.81) converges to the solution of the 
original Poisson equation ed, i- e - Theorem 11.11 in Section fTTTl To prove this theorem, we split it to two 
parts. First, we prove that the solution of the integral equation (11.61) converges to the solution of the Poisson 
equation ed, which is given in Theorem 13.21 Then we prove Theorem l3.3l to show that the solution of (11.81) 
converges to the solution of (11.61) . 


3.1 Integral approximation of Poisson equation 

To prove the convergence of the integral equation (11.61) . we need following theorem about the consistency 
which is proved in I20j . 

Theorem 3.1. Let u(x) be the solution of the problem ed- Let u £ H 3 ( At) and 

r(x) = j f R t (x, y)(u(x) - u(y))dy - [ i? t (x, y)/(y)dy. 

1 J m Jm 

There exists constants C, To depending only on M and dM, so that for any t <Tq, 

IK x )IIl 2 (a ^) < Ct 1/2 \\u\\ H3{M) , (3.1) 

l|Vr(x)|| i2(M , ) < C\\u\\ H3{M) . (3.2) 


Using the consistency result, Theorem 13.11 and the stability results presented in Section [2j we can get 
following theorem which shows the convergence of the integral equation (11.61) . 

Theorem 3.2. Letu(x) be solution of ED and Ut(x) be solution of ED- There exists C > 0 only depends 
on At and dM, such that 


ll w - u t \\m{M' t ) < ct^wn m(M) 




Proof. Let et(x) = w(x) — u t (x), first of all, we have 

7 [ e t (x) [ R t (x, y)(e t (x) - e t (y))dydx (3.3) 

1 J M' t J M 

= 7 / e t (x) / i? t (x,y)(e t (x)-e t (y))dydx+i / e t (x) / i? t (x,y)(e t (x) - e t (y))dydx 
1 JM' t JM' t 1 JM't JV t 

[ [ R t (x, y)(e t (x)-e t (y)) 2 dxdy+i / e t (x) [ i? t (x, y)(e t (x) - e t (y))dydx. 

jm: Jm' 1 J Mi jv t 


JM’t Jm 

The second term can be calculated as 
1 


f e*(x) f R t (x, y)(e t (x)-e t (y))dydx (3.4) 

^ JM' t JVt 

7 / |et(x)| 2 ( [ R t (x, y)dy N ) dx- j f e t (x) ( f R t (x,y)u(y)dy 

1 JM' \JVt J 1 JM' \JVt J 


Here we use the definition of e* and the volume constraint condition u*(x) = 0, x G Vt to get that e*(x) = 
u(x), x e Vt. 

The first term is positive which is good for us. We only need to bound the second term of (13.411 to show 
that it can be controlled by the first term. First, the second term can be bounded as following 


< 


< 


< 


< 


< 


[ e *( x ) ( [ i? t (x,y)w(y)dy N ) dx 

JM't \Jv t J 

\ Jm, M x )I (f v - R t(x,y)dy^j ^ i? t (x, y)|w(y)| 2 dy^ 

\ 7l e t( x )| 2 (^j v ^t(x,y)dy^ dx + 2 ^ i? t (x, y)|u(y)| 2 dy 

Y t J m , l e *( x )| 2 ^(x,y)dy) d x+J J v l M (y )| 2 #t( x , y ) dx j dy 

Y t J M , l e *( x )| 2 ^(x,y)dy) dx+fjf Ky)| 2 dy 

^ Jm' ' e ‘^ 2 (X i?t( ' x,y ' )dy ) -^v.ll/ll 2 !^). 


(3.5) 


Here we use Lemma lA.il in Appendix A to get the last inequality. 
By substituting (13.51) . (13.41) in (13.31) . we get 


> 


/ e t (x) / tft(x,y)(e t (x)-e t (y))dydx 
JA4J Jm 

77 / [ Rt{x, y)(e t (x) - e t (y)) 2 dxdy 

I M ' (X i?t< ' x,y ^ dy ) dx_c 'ii/n^ i (x)^- 


(3.6) 


1 

_ 2f 


This is the key estimate we used to get convergence. 
Notice that e*(x) satisfying an integral equation, 


[ #t(x,y)(e*(x) - e*(y))dy = r(x), VxeM' ( , 
jm 


(3.7) 
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where r(x) = \ f M R t (x, y)(u(x) - u{ y))dy - f M Rt(x, y)/(y)dy. 

From Theorem 13.II we know that 

ll r ( x )ll,L 2 (At;) - 

|| Vr (x) || ^ 2 (^ /) < C\\u\\ H 3( M ) < C\\f\\ HHM) . 


Now we can get L 2 estimate of e*. Using Corollarv l2.2l we have 


W e t\\lM’ t -T f f • R t( x -y)( e t( x )-e t (y)) 2 dxdy 
1 JM' t JM' t 

+ j j , l e *( x )| 2 (yj -Rt(x,y)dy^ dx 

/ e t (x) / i?t(x,y)(e t (x) - e t (y))dydx 
J Mi J M 


(3.8) 

(3.9) 


(3.10) 


(from (13.61) 1 < 


C 

7 




(from E3) < C||et|| 2 ) m;|| r|| 2 ,Mj + C\\f\\ 2 m{M) Vt 

(from j3U)) < C||/|| ff i (A 4 ) ||e t || 2 ,^;V / t + C'||/||^i ( _ M) v / t. 


This gives that 


lletlb.Ai; < C't 1//4 ||/|| ff i( 2 v ( ). 


(3.11) 


Next, we turn to estimate the L 2 norm of Ve* in A4' t . Using the integral equation (13.71) . et has following 
expression 


If t 

et(x) = —-- / 7?t(x, y)e t (y)dy H ——r(x) 

wt{x) w U x ) 

[ Rt(x, y)e t (y)dyH- [ R t {x,y)u(y)dy H-^-r(x). 

Cvti w H x ) dV t W H X ) 


(3.12) 


Wt(x) 


Then ||N7e t ||| can be bounded as following 


IVetHl^j <C 
1 


+U 


w t (x) 


'v t 


Wt (x) 
i? t (x,y)w(y)dy ) 


f i?t(x,y)e t (y)dy) 

At; / 


(3.13) 


2 ,At' 


Ct 2 


2 ,At; 


K x ) 

W t (x) 


2 .At; 


Corollary 12.11 gives a bound the first term of (13.131) . 

2 

/i r \ 

V 


—rw [ Rt (x, y)et(y)dy j 
w t (x) J M , t J 


(3.14) 


2 , At; 


< 7 -/ [ Rt(x,y)(e t (x) - e t (y)) 2 dydx+ y [ |e t (x)| 2 ([ R t {x, y)dy 

1 JM' t JM' t 1 JM' t \Jv t 

The second and third terms of (13.131) can be bounded by direct calculation. 


dx. 


v( r( *() 

2 

< c 

Vr(x) 

2 

+ c 

r(x)Vro t (x) 

\wt(x)J 

2,M't 

w t (x) 

2 ,At; 

M x )) 2 


2 ,At; 


(3.15) 


C 


< C ||Vr(x)|| 2i _ M , + — ||r(x )|| 2 _ M , 

< C|l/ll?fi(At)’ 
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and 


< 


w t (x) J Vt 
Vfflt(x) 


K(x )) 2 


IVt 


i? t (x,y)w(y)dy ) 
i?t(x,y)it(y)dy 


+ 


w t (x) 


'v t 


Vi?t(x, y)w(y)dy 


< C'||/|| ff i ( 7 K ) f R t (x,y)dy + C||/|| ff i (7W) [ |i?' t (x,y)|dy. 

Jv t Jv t 


Then the second term of (13.141) has following bound 


i?t(x, y)w(y)dy 


wt(x) Jv t 


i.,M' 


< G \\f\\ 2 HHM) J m/ (J v R t(x,y)dy S j dx + C\\f\\ 2 H1{M) ^ |i?t(x,y)|dy 

< c||/lllfi(M) / / ^t( x .y) d y dx + c 'll/ll^i(>i) / / ^(x.y^dydx 

jaij Jv t Jm; Jv t 

Now we have the bound of ||Ve t || 2 .A^; by combining (13.131) . (13.151) . (13.141) and (13.171) 

7 -I [ R t(x, y)(e t ( x ) - e t (y)) 2 dxdy 

1 jjm; ja'i; 

|et ( x )| 2 ( [ R t (x,y)dy\ dx + C\\f\\ 2 m{M) Vt. 
i' t \JVt J 


dx 


|Ve t ||i M , < 


t JM' 


Then the bound of ||Ve t || 2 ,x; can be obtained also from (13.181) 

C 


\\Ve t r 2M , <- 


(from (13.61) ) < 


[ [ Rt(x,y)(e t (x) - e t (y)) 2 dxdy 

JM' t JM' t 

b J J M , l e ‘( x )| 2 (f v R t(x,y)dy S j dx + C\\f\\ 2 H1{M) Vi 

7 - / e t (x) f Rt(x,y)(e t (x) - e t (y))dydx + C'||/||hi ( x ) 
1 J M' J M 


Vt 


(from (EH)) <\\eth,M' t \\ r h,M' t + C\\f\\ 2 H i {M) Vt 
(from (EU)) <C\\f\\ H i {M) \\e t \\ 2 , M[ Vi + C\\f\\ 2 Hl{M) Vi 

(from (13.11|)) <Ct 3 / 4 \\f\\ H i (M) + C\\f\\ 2 m(M) Vt. 


Then we have 


l|Vet|| 2) x; < C't 1 / 4 ||/|| ff i( A1 ). 

The proof is completed by putting (13.111) and (13.201) together. 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 

□ 


3.2 Discretization of the integral equation 

Suppose u = [ui, • • • , UnY is the discrete solution which means that it solves (11.81) . First, we interpolate the 
discrete solution from the point cloud P = {pi, • • • , p„} to the whole manifold M. Fortunately, the discrete 


11 




































equation (11.81) gives a natural interpolation. 


u t , h (x) = < wt,h{ x ) 


Y Rt(x,Pj)ujVj +t Y R t{*,Pj)f(Pj)Vj ) » x EK 

p j&M 

0 , 


(3.21) 


x e v t . 


Then, we have following theorem regarding the convergence from ut,h to ut- 


Theorem 3.3. Let u t (x) be the solution of the problem m and, u be the solution of the problem la. if 
f € C 1 (M ) in both problems, then there exists constants C > 0 depending only on Ai and dAi so that 

\\u t ,h ~ ut\\m(M' t ) < ^■||/llc 1 (M)i 

as long as t and are both small enough. 

To prove this theorem, we need the stability result, Theorem 12.11 and the consistency result which is 
given in Theorem 13.41 

To simplify the notations, we introduce some operators here. 


L t u{x) = \ ( Rtf*-, y)( u ( x ) — w(y))dy + —^ f i? t (x,y)dy, 

1 JM’ t 1 Jv t 

L' t u(yf) = - [ i?t(x,y)(u(x)-u(y))dy, 

1 JM' t 

2 ^ u(x.) r 

{L t , h u)(x) = - Y R t(x,Pj)(u( x ) ~u(Pj))Vj + —— Y R t(Pi’Pj) V L 


P j£M' t 


Pj&Vt 


( L 't,h u )(*) = 7 E Rt(x,Pj)(u(x.) - u(pj))Vj. 


Pj £M' t 


(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 


It is easy to check that u t ,h satisfies following equation if u = [ui, ■ ■ • , u n ]* solves (11.81) . 

Z! R t(^Pj)f(P,)Vj 

PjGM 

Now, we can state the consistency result as following. 

Theorem 3.4. Let iq(x) be the solution of the problem (toil and u be the solution of the problem da- u 
f € , in both problems, then there exists constants C > 0 depending only on Ai and dAi so that 


Ch 


\\L t {u t ,h ~ ut) IIl 2 04 ) < fi/zWfWc^M) 

Ch 

\^Lt{ut,h- U t )\\L 2 (M' t ) < -p-WfWc 1 ^)- 


(3.27) 

(3.28) 


as long as t and -^= are small enough. 


4 Convergence of the eigenvalue problem 

In this section, we investigate the convergence of the eigenvalue problem (11.111) to the eigenvalue prob¬ 
lem (11.91) . First, we introduce some operators. 

Denote the operator T : L 2 (M) —► H 2 (Ai) to be the solution operator of the following problem 

f A M u(x) = /(x), x e M, 

\ u(x) = 0, x e dAi. 
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where n is the out normal vector of A4. 

Denote T t : L 2 {M) L 2 {M) to be the solution operator of the following problem 

f -7 / #t(x,y)(u(x) - u(y))dy = [ y)/(y)dy, x € M' t 

N 1 JM JM 

[ u(x) =0, x G Vt 

The last solution operator is T tt h ■ C(M) —> C{Ai) which is defined as follows. 


Tt,h(f)( x ) = < w t>h {-x) 



i? t (x, p j)UjVj - t 


PjGM 


X G M' 


y 0, x g v 4 , 

where iu tj /,(x) = Y^ Pj & m ^t( x i P j)Vj and u = (iti, • • • , u n )* with Uj = 0, p_, G Vt solves the following linear 
system 


1 

t 


Y Rt(Pi,Pj)(ui 

PjGM 


~ u j) v i= Y R t(Pi’Pj)f(Pi) V r 

Pj£M 


We know that T, T t and T ti h have following properties. 

Proposition 4.1. For any t > 0, h > 0, 

1. T, T t are compact operators on H 1 { At) into H 1 { At); T tl T t ,h are compact operators on C 1 (M) into 

C X {M). 

2. All eigenvalues of T,T t ,T t ,h are real numbers. All generalized eigenvectors of T,T tl T tt h are eigenvec¬ 
tors. 


Proof. The proof of (1) is straightforward. First, it is well known that T is compact operator. T t) h is actually 
finite dimensional operator, so it is also compact. To show the compactness of T t , we need the following 
formula, 

T t u =—r^l R t(*,y) T My)<iy H-fw [ R t {x.,y)u(y)dy, Mu&H 1 {M). 

«M X ) Jm w H x ) Jm 

Using the assumption that R G C 2 , direct calculation would gives that that T t u G C 2 . This would imply 
the compactness of Tt both in H 1 and C 1 . 

For the operator T, the conclusion (2) is well known. The proof of T t and T t ,h are very similar, so here 
we only present the proof for T t . 

Let A be an eigenvalue of T t and u is corresponding eigenfunction, then 


L t T t u = XL t u 


which implies that 

A = f M f M -Rt(x,y)u*(x)u(y)dxdy 

f M u *( x )( L t u )( x '> dx 

where u* is the complex conjugate of u. 

Using the symmetry of L t and R{x., y), it is easy to show that AgI. 

Let u be a generalized eigenfunction of T t with multiplicity m > 1 associate with eigenvalue A. Let 
v = (T t — A ) m ~ 1 u, w = (T t — A ) m ~ 2 u, then v is an eigenfunction of T t and 

T t v = Xv, {T t — X)w = v 
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and v(x) = 0, w(x) = 0, x € V t . 

By applying L t on both sides of above two equations, we have 

A L t v = L t (T t v) = [ Rt(x, y)u(y)dy = [ R t (x,y)v(y)dy, x e M' t 
J m J M' t 

L t v = L t (T t w) — XL t w = / R t (x,y)w(y)dy - XL t w, x e M' t 

J M' t 

Using above two equations and the fact that L t is symmetric, we get 

0 = ( w, XL t v — / R t (x,y)v(y)dy 

\ 

= (XL t w- / R t {x,y)w{y)dy,v 

\ J M' t 

= i L tv,v) M ' t > C\\v\\l 

(T t — X) m ~ 1 u = v = 0. This proves that u is a generalized eigenfunction of T t with 
Repeating above argument, we can show that u is actually an eigenfunction of T t . □ 

Theorem 13.21 actually gives that T t converges to T in H 1 norm. 

Theorem 4.1. Under the assumptions in Section \l.R for t small enough, there exists a constant C > 0 
such that 


which implies that 
multiplicity m — 1. 




\\T-T t \\ m < Ct x '\ 

Using the arguments in [20] and Theorem 13.81 we can get that T t: h converges to T t in C 1 norm. 

Theorem 4.2. Under the assumptions in Section Q77J for t , h small enough, there exists a constant C > 0 
such that 

C!h 

\\Tt,h - T t \\ C i < fd/I+2' 

And we also have the bound of T t and T f j L following the arguments in E01 . 

Theorem 4.3. Under the assumptions in Section li.il for t , h small enough, there exists a constant C 
independent on t and h, such that 

\\Tt\\m < C, Halloo < Ct~ d !\ ||T t , h || C i < 

. Before state the main theorem of the spectral convergence, we need to introduce some notations. Let 
A be a complex Banach space and L : X —> X be a compact linear operator. p(L) is the resolvent set of L 
which is given by z £ C such that z — L is bijective. The spectrum of L is cr(L) = C \p(L). If A is a nonzero 
eigenvalue of L, the ascent multiplicity a of X—L is the smallest integer such that ker(A— L) a = ker(A— L) a+1 . 

Given a closed smooth curve T C p(L) which encloses the eigenvalue A and no other elements of <r(L), 
the Riesz spectral projection associated with A is defined by 

£(A,L)= 2 h ~ L) ~ ldZ ' 

where i = y/—l is the unit imaginary 

Now we are ready to state the main theorem about the convergence of the eigenvalue problem. And its 
proof can be given from Theorems 14.1114.21 and 14.31 following same arguments as those in J2D3 • 
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Theorem 4.4. Under the assumptions in Section \l.ll let A, be the ith smallest eigenvalue of T counting 
multiplicity, and be the ith smallest eigenvalue ofT t counting multiplicity, then there exists a constant 
C such that 

I Xt i h ’ 

and there exist another constant C such that, for any (f> € E(Xi,T)X and X = H l (AA), 

U-E{al' h ,T t M m(M) < C (V/ 4 + . 

where a\' h = {A‘ ,/l e cr(T tth ) : j € h} and I, = {j £ N : A j = AJ. 

The convergence result, Thorem ll.2l follows easily from the above theorem and Proposition 14.II 


5 Numerical results 

In this section, we present several numerical results to show the convergence of the Point Integral method 
with volume constraint, PEVLVC for short, from point clouds. 

The numerical experiments were carried out in unit disk. We discretize unit disk with 684, 2610, 10191 and 
40269 points respectively and check the convergence of the point integral method with volume constraint.In 
the experiments, the volume weight vector V is estimated using the method proposed in m- First, we 
locally approximate the tangent space at each point and then project the nearby points onto the tangent 
space over which a Delaunay triangulation is computed in the tangent space. The volume weight is estimated 
as the volume of the Voronoi cell of that point. 

Table[l]gives the l 2 error of different methods with 684, 2610, 10191 and 40269 points. The exact solution 
is cos 2tt^/x 2 + y 2 . PIM_Robin is the Point Integral method and using Robin boundary to approximate the 
Dirichlet boundary condition, i.e. solving the integral equation (11.511 and here /3 is chosen to be 10 -4 . PIM_VC 
is the Point Integral method and using volume constraint to enforce the Dirichlet boundary condition. These 
two methods both converge. The rates of convergence are very close and the error of PIM_VC is a little 
larger than the error of PEVLRobin. 


1^1 

684 

2610 

10191 

40269 

PIM_Robin 

0.1500 

0.0428 

0.0140 

0.0052 

PEVLVC 

0.3046 

0.0747 

0.0201 

0.0067 


Table 1: l 2 error with different number of points. FEM: Finite Element method; PEVLRobin: Point Integral 
method with Robin boundary; PEVLVC: Point Integral method with volume constraint. The exact solution 
is cos27r-\/a; 2 + y 2 . 

Fig. [2] shows the result of the eigenvalues of Laplace-Beltrami operator with Dirichlet boundary in unit 
disk. Clearly, the eigenvalues also converge and the larger eigenvalues have larger errors which verify the 
theoretical result, Theorem 11.21 

Above numerical results in unit disk are just toy examples to demonstrate the convergence of the Point 
Integral method with volume constraint. However, our method applies in any point clouds which sample 
smooth manifolds. Fig. [3]shows the first two eigenfunctions on two complicated surfaces (left hand and head 
of Max Plank). 

6 Conclusion 

In this paper, we use the volume constraint [7] in the Point Integral method to handle the Dirichlet boundary 
condition. And the convergence is proved both for Poisson equation and eigen problem of Laplace-Beltrami 
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Figure 2: Eigenvalue given by Point Integral method with volume constraint in unit disk. 


operator from point cloud. Our study shows that Point Integral method together with the volume constraint 
gives an efficient numerical approach to solve the Poisson equation with Dirichlet boundary on point cloud. 
In this paper, we focus on the Poisson equation. For other PDEs, we can also use the idea of volume 
constraint to enforce the Dirichlet boundary condition. The progress will be reported in our subsequent 
papers. 

Acknowledgments. This research was partially supported by NSFC Grant (11201257 to Z.S., 11371220 
to Z.S. and J.S. and 11271011 to J.S.), and National Basic Research Program of China (973 Program 
2012CB825500 to J.S.). 


A One basic estimates 


Lemma A.l. Let u(x) be the solution of m and f € H 1 (Ai), then there is a generic constant C > 0 and 
To > 0 only depend on Ai and dAi, for any t < To, 


>v t 


Ky)| 2 dy < Ct 3 / 2 \\f\\ 2 Hl(M) . 


Proof. Both Ai and dAi are compact and C°° smooth. Consequently, it is well known that both Ad and 
dAi have positive reaches, which means that there exists To > 0 only depends on At and dAi, if t < T 0 , V t 
can be parametrized as (z(y), r) £ <9Af x [0,1], where y = z(y)+r(z'(y)— z(y)) and det ^ d ( z ^ <C\ft 
and C > 0 is a constant only depends on Ai and dAi. Here z'(y) is the intersection point between dAi' 
and the line determined by z(y) and y. The parametrization is illustrated in Fig|4] 
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Figure 3: The first (upper row) and second (lower row) eigenfunctions with Dirichlet boundary. 



First, we have 


dy 


[ |u(y)| 2 dy= [ |u(y) - u(z(y))| 2 dy 
Jv t Jv t 

ff 1 d 2 

= Jv J ds u ^ y + s ( z ( y ) - y )) ds d y 

= [ I (z(y) -y) ■ Vw(y + s(z(y) - y))ds 
Jv t Jo 

< ct [ f |Vu(y + s(z(y)-y))| 2 dsdy 

Jv t Jo 

< Ct sup [ |Vu(y + s(z(y) — y))| 2 dy. 

0 <s<l Jv t 

Here, we use the fact that ||z(y) — y ||2 < 2 yft to get the second last inequality. 
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Then, the proof can be completed by following estimation. 


/ |Vit(y + s(z(y) — y))| 2 dy 
Jv t 

< CVt [ [ |Vw(z(y) + (1 - s)r(z'(y) - z(y)))| 2 dz(y)dr 

JO JdM 

< CVt sup [ |V?i(z + (1 — s)t(z' — z ))| 2 dz 

0<r<l JdM 

< CVt sup f |Vu(z)| 2 dz 

0<r<l Jr BiT 

— CVt\\ u \\ 2 H -2(M) < 


where T SjT is a k — 1 dimensinal maniflod given by T SjT = {z + (1 — s)r(z' — z) : z £ ctAd}. We use the trace 
theorem to get the second last inequality and the last inequality is due to that u is the solution of the Poisson 
equation (fTTH . □ 


B Proof of Theorem 13.4 

First, we need following important lemma which tells us that the discretized scheme is stable in l 2 sense. 

Lemma B.l. For any u = ,u n Y with Ui = 0, p, £ Vt, there exist constants C > 0, Co > 0 

independent on t so that for sufficient small t and 

J J2 RtiPiiPjKui-UjfViVj u i V f 

PiGM Pj£A4 v pjGM 


Proof. First, we introduce a smooth function u that approximates u at the samples P. 

1 


(x) = -— Y R t' ( x , P*)uiVi, x £ M, 

wt',h(x) frf 


(B.l) 


where Wt\h{ x ) = and t' = t/ 18. Using the condition that Ui = 0, Pi £ V t and 

t' = t/18, we know that 


it(x) = 0, x £ Vt'. 

Then using Lemma 12.21 we have 

f |w(x)| 2 dx < — f I R t ' (x, y) (m(x) — u(y)) 2 d/r x d/r- 
J M 1 J M J M 


(B.2) 
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On the other hand 


■'52Rt>(x-,Pi)uiVi - —— ^ Rt> (p, , y)u 3 V :l | d/i x d/i y 


w t'A x ) ’ w t'Ay)~[ 


/ / Rf (x, y) ('u(x) - u(y)) 2 d/i x d/i y 

JA1 

= / / ^t'(x,y) 

JM JM 

= LL R, ' {x ’ y) (»,, k (*K- A (y) E i ^ ( X ,P«)A'(P i ,y ) KVi ( «,-»,) J ta*, 

<[ [ R t ’ (x, y)---i-- V i?t/(x,p l )i? t /(p i ,y)^y j (M i - u^^^d^y 

JmJm W t ',h{X-)W t > ,h{y) j~x 

= Y] ( I [ -TT- 

^ 1 \JMJM w t',h(K)Wt’,h{y) J 


Denote 


A = 


1 


-i? t / (x, pO-Rt' (Pj, y)i?t' (x, y)d/x x d/iy 


/A4 /m wt>,h(x)w t >, h { y) 

and then notice only when |pj — pj | 2 < 36£' is A ^ 0. For |p; — pj | 2 < 36£', we have 

| 2 \ - 1 


A < 


< 


M J M 

CCt 


R t ' (x, p i)R t ' (p j, y)R t ' (x, y)R 


I Pi -Pjt 


JO JM JM 


Rf ( x , Pi)R t ' (pj,y)R 


<CC t / / R t '(x 1 Pi)R t '{P],y)R 

JM JM 


72 t' 

|Pi ~ Pjl 2 

72 f 

P» - Pjl 2 

72 t' 


R 


I Pi - PiY 


72 t' 


d/i x d/iy 


d/J X d/iy 


dfi x dfi y < CC t R 


(^) 


Combining Equation (IB.31) . (IB.41) and Lemma 12. 2 1 we obtain 

V Rt (PuPj) {Ui - Uj) 2 ViVj > / |u(x)| 2 d/r x 

p,- — Jm 

Denote 


Pi.p jSM 


B = 


C t 


R 


|x-Pi 


'm 


n ("i 

E 


™l, h (Pj) V 4t' 


At' 

Ipj - Pi 1 2 


I*- Pip 

At' 


d/i x 


|Pj -Pi\ : 

At' 


Vj 


(B.3) 


(B.4) 


(B.5) 


and then \B\ < jtji- At the same time, notice that only when |p* — P /| 2 < 16£' is B ^ 0. Thus we have 

|B| 
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and 


r. n 

/ M 2 (x) d/! x ~' Y ' u 2 { pj ) V j 
J M j=1 


(B. 6 ) 


< ^ \C t UiUiViVi\\A\ 

i,l— 1 


c 

- + 1/2 


i,l— 1 


C t i? 


|p* - p i\ : 

32 t' 


UiUiViVi 


^ Ch r , D f | Pi — Pz | 2 ^ 2 T /T/ ^ ^ 2 T / 

- t l/2Z^ Cti? l 39*/ ) UiViVi< +1/2 2^«iVi. 


2 ,Z =1 


32f' 


t !/2 / 


i=l 


Now combining Equation (IB. 51) and (IB.61) . we have for small t 


^2u 2 {p t )Vi = 


< 


/ M 

JM 

CCt 


Ch _ " . 

(x)d/X x + Tin J2 U i V i 


t l/2 

t, R <“< - 

2,j = l ' ' 2=1 


2 — 1 
2 


£ 


Here we use the fact that for t = ISt' 


R 


Pi-Pj| 2 \ ^ 1 T,f\Pi-Pj\ 


At' 


< t r 

Oo 


At 


Let 5 = 2 u, If ELi u 2 (Pi)^ > 52 EL i u 2 Vl, then we have completed the proof. Otherwise, we 


have 


y^Xuj - u{pi)) 2 Vi = ^u 2 Vi + ^u(p i ) 2 V i - 2^ruiu{pi)Vi > (1 - 

2=1 2=1 2=1 2 = 1 2=1 

This enables us to prove the cllipticity of C in the case of Y^i= i u 2 (P*)^* < ELi u i ^ as follows. 


I P» ~ Pjf 

At' 

Ip* - Pjf 

At' 


{m - UjfViVj 


Ui(ui - Uj)ViVj 


c t j2 R 

*0 = 1 
n 

= 2 C t J2 R 

*o=i 

n 

= 2 ^uiiui - u(pi))w t ,h(Pi)Vi 
2=1 

n n 

= 2 ^2(ui - u(pi)) 2 w t ,h(Pi)Vi + 2^u(p;)(uj - u(pi))w t ,h(Pi)Vi 

2=1 2=1 

n / n \ / n \ 

> 2- u(Pi)) 2 w t ,/ l (Pi)^ / * ^ 2 ( ^u 2 (p i )w ti / l (p i )y i ) ( ^(ui - u(Pi)) 2 w ti h(p i )’l / i ] 


2=1 


\i=l 


\ 2 = 1 


n / n \ / n \ 

zL 2n7 m in ^ ^ (^2 ^(Pi)) 2u; max | ^ ^ U (p i)Vi J | ^ ^ (^2 ^(Pi)) ^ ] 


2=1 


vi=l 


V2=l 


> 2(io min (l - S) 2 - w max S( 1 - (5)) ^ u 2 Ei > w min (l - S) 2 ^2 u i v i- 


2=1 


2=1 


□ 
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One direct corollary of above lemma is the boundness of (E Pi e.M u i Vi) 


1/2 


Corollary B.l. Suppose u = (ui, ■ ■ ■ ,u n y with m = 0, p; G Vt solves the problem (11.81) with f G C(M). 
Then there exists a constant C > 0 such that 

\ 1/2 

E u i v *) ^ 

K PiGM 

provided, t and are small enough. 

Proof. From the elliptic property of C, we have 


E u ^ Vi -EE uiVi 


1/2 


/ 

E 

ypieM 


> E R t(Pi’Pj) V i ) V i 

PjGM 


1/2 


PiGM piGM \PjGM 

< (E u ? Vi 

\PiGM 

< c ( e 

\Pi GM 

This proves the lemma. 

We can also get the bound of ut,h which is defined in (13.211) as following 


1/2 


□ 


1 


Ut,h(x) = < Wt,h(x ) 


E R t(x,Pj)ujVj +t E R t(.X’Pj)f(Pj) V j\ » xeM't, 


i PjGM 


PjGM 


0 , 


x G V t . 


Lemma B. 2. Let u= ,«„]* be the solution of the problem (11.81) with f G C(A4) andut^h be associate 

smooth function defined in (13.211) . Then there exists C > 0 such that 


\\ut,h{x)\\ L 2 {M , t) < CH/lloc 

l|V«t, fc (x)|| i2(A<i) < ^ll/ll 
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Proof. First, ||u t ,/ l (x)|| i2(A<( , ) is bounded 

||'Wt,7i( x )|| i 2(JV(/) 

f 1 


lM' t ™£h( X ) 


Y R t( x ’Pj) u jVj-t Y R t(*,Pj)fjVj ] dx 

Pi 6 p PiGP 

2 

2 


< C I I Y2 R t{ x iPj) u jVj ] dx + ct 

yPj£P 




>M' 


Y R t( x ,Pj)fjVj ] dx 

Pi 6 p 


; c 




YY Y2 R t fx.,Pj)u?Vj J dx + Ct 


2 ii ni 2 
I oo 


, PjGP 


iPjEP 


IM' 


Y R t( x -pi) v J ) dx 

yPjGP 


< C Y u 2 V 3 / P t (x, Pi )dx +Ct 2 

P, e p V M 't J 

< cY 2 u2 i V i + Ct2 \\f\\lo < CWfWlo 

Pi 6 P 


Using similar arguments, we can get the bound of || ’Vu tt h(x-)\\ L2 ( M ,y From the definition of Utji- we can see 
that all derivatives are applied on the kernel functions. The kernel functions are smooth functions, it gives 
one factor of after derivative. □ 

Now, we are ready to prove Theorem 13.41 


(B.7) 


Proof, of Theorem \3.4\ 

First, we split L t (ut,h ~ ut) to three terms, for any x £ A4' t , 


L t ( u t ,h ~ u t ) 

— R t(jl‘t,h') R t.h{ut.h) T I J t,h{.'U‘t,h) R tV^t 

U t ,h{ x ) 


— Yt( u t,h) R t,h( U t,h)) + 


>Vt 


R t( x , y)dy — Y R t( x ,Pi) v i 


PiGVt 


+ ( / 7? t (x, y)/(y) — Y R t( x iPj)f(Pj)Vj 


pjGM 


To get the last equality, we use that ut and Ut } h solve equation (11.Gl) and equation (13.261) respectively. 

The second and third terms are easy to bound. By using Lemma IB.21 and (P, V) is h-integrable approx¬ 
imation of A4, we have 


U t ,h( x ) 


>v t 


R t( x , y)dy- Y Pt(x, P,;)U 


PiGVt 


U t ,h( x ) 


/V t 


Pt(x,y)dy- Y Pt(x,Pt)U 


Pi eVt 


Li(M’t) 


L 2 (M' t ) 


< 


< 


£3/2 11^11°°’ 

(B.8) 

Ch . 

-pWfWoo 

(B.9) 
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and 


[ #t(x,y)/(y)- Y i? t( x >Pj)/(Pi) 1/ i 

JM Pj gm 

R t( x ,y)f(y)~ Y R t( x ,Pj)f(Pj) v i 

P i&M 


L 2 (M' t ) 


Ch 

< ^WfWcHM), (B.10) 


) M 


Ch 

< —WfWcHM) (B.ll) 


L 2 (M‘ t ) 

The first term of IB.7)1 is much more complicated to bound. We split it further to two terms. Denote 

1 


a t ,h{x) = 


w t ,h { x ) 


Ct,h ( x ) = - R t( x -Pj)f(.Pj) V ji 

w tM x ) 


Y R t( x -P:i) u 3 V v, 


PjGM' t 


PjGM' t 


and then = a t ,h{ x ) + ^^(x), x e .AdJ. 

First we upper bound \\L' t (u t , h ) ~ L' t h {u t ,h)\\L 2 (M' t )- For °t,h, we have 


| ( F i c t,/i R t,h c t,h) ( x ) | 

/ i?t(x,y)(ct, ft (x) - c t , h (y))dny - V i?t(x,pj)(c t ,h(x) - c t ,h(Pj))Vj 
JM' 


PjGM't 


< -\ct,h(x)\ 


IM 


i? t (x,y)d/i y - Y R t( x - P:i) V i 

Pj€M' t 


+ - 


Ch 


/ i? t (x, y)c t ,ft(y)d/x y - V R t (x, Pj)ct, h {Pj)Vj 

JM' rt.CZ 


Pj£M' t 


Ch 


— j.3/2 l C *. ?1 ( X )l d” ^3/2 11 IIC 1 (7W') 

Ch . Ch 

*372 




Ch. 


+ £ 1 // 2 ||/I|oo) < —IU lloo- 


For at,h, we have 


(at./iW) 5 


< 
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J M' 

Ch 2 
t 

Ch 2 
t 

ch 2 r 

t J M 

Ch 2 


/At' 


-Rt( x , y)d/i y — ^ /C(x,p 7 )V,- 

PjGAIJ 


d// x 


(a t ,h( x )) d/x x 


'At' 


'At' 


w t ,h(y 


Y Rt( x ,Pj)ujVj dp x 
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Rtfr’Pj) u j v 3 R t( x ’Pi) V o d/i x 


VP) 6 p 
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Ch 2 


< — | ^ w 2 ^- I i? t (x, Pj )d/z x j <— ^ 


. Pj GAtJ 


/At 
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(B.14) 
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Let 


A = C t 

-<* E 
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m’ w t,h(y) 
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x-y r 
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Pj&M 


, w t ,h(Pj) \ At 


4 1 

|x - P j 


I Pi - yr 


4 1 


d^iy 


R 


I Pi - Pjl 

At 


Vi- 


We have \A\ < for some constant C independent of t. In addition, notice that only when |x — p;| 2 < 16i 
is A ^ 0, which implies 


\A\<-\A\R 

oo 


|x~Pd 
32 1 


Then we have 


/ / R t (x.,y)a t , h (y)diJ,y - E R t(^Pj)a t ,h{P]) v j 

Jm ' Jm ' Pj eM' t 

E CtUiViA ] dfi 
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<- 
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ctf_ 

t 

,ctf_ 

t 
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\Pi&M' t 
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'M' 


E Ct\ui\ViR 

\Pi&M' t 


|x-p,| 

32t 


E C ‘ R 

\P i£M' t 


|x - P,: 
32 1 


nfVi 


d/i x 


E C < R 


|x - P, 
32t 


bi I d/i x 


E u ^ Vi 




C t R 


’ M' 


|X - P, 
32 1 


d fi x < 


Ch 2 


E u i Vi 

\Pi£M' t 


Combining Equation (IB. 141) . (IB. 151) and Theorem lB.il 

II L’t a t,h — L' th a ti h\\L 2 (M) 

= If \{Lt(a t< h) ~ L' th (at,h)) (x)| 2 d/i x 
V JM' 


1/2 


1/2 


Ch 


Ch 


- £3/2 ( E ) - *3/2 


£3/2' 


yPieMJ 

Assembling the parts together, we have the following upper bound. 

II L t Ut,h ~ ^t,h u t,h\\L 2 (M' t ) 

— I|-^'t a t,ii ^t,h a t,h || L 2 (M' t ) T ll^'t c t,ii ^t,h c t,h || L 2 (M' t ) 


< 


Ch 
*3/2 I 


+ 


Ch, 


Ch | 

l°° — £3/2 I 


(B.15) 


(B.16) 


The complete L 2 estimate follows from Equation (IB.8I) . (IB.1QI) and (IB. 161) . 

Next, we turn to upper bound ||V(L' t u t - L' t h u t ,h)\\L 2 (M' t )- Consider \\S7(L' t a t ,h - L' t h a t,h)\\L 2 {M' t ), it 
can be splited into the summation of three terms. Next, we estimate these three terms separately. The first 
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term is 
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(B.17) 


where i? 2 t(x, p^) = C t R ^. Here we use the assumption that i?(s) > do for all 0 < s < 1 / 2 . 

The second term is 


< 


Ch 2 

? 2 


Iat.^x)! 2 / Vi? t (x,y)d/r y - E Vi? t (x, p^V,- 
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/ |< 2 t,/i(x)| 2 clp x < — 55 
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Let 
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I x - Pj 
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I Pi - yf 
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Pi ~ Pj 
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d/iy 
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We have |B| < ^75 for some constant C independent of t. In addition, notice that only when |x — x ^| 2 < 16? 
is B ^ 0, which implies 
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Then we have the upper bound of the third term 


/ / VJ? t (x,y)a t ,/ l (y)d/x y - V Vi? t (x,p J )a tih (p J )V 7 

Jm' Jm' 
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Combining Equation (IB. 171) . (IB. 181) and (IB.191) . we have 

II V(L t at,ft, — L th at,h)\\L 2 (M' t ) 


1/2 




| (L t {a tt h) - L tih (a t ,h)) (x)| 2 d/r x 


1/2 
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Using a similar argument, we obtain 

\\^(Lt c t,h — -^t,h c t,h)\\L 2 (M' t ) — £3/2 

and thus 

II N(L t Ut,h — L th Ut,h)\\L 2 (M' t ) — ~j% 

At last, we complete the proof using (IB. 91) . (IB.Ill) and (IB.201) 


Ch, 


d/i x 


(B.19) 
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